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1. Abstract 

T Se tank Str °" gly affcct the “ nditi » of*, tank 
ana its ability to withstand operational stresses. These thermal stresses also affect the 

growth of any surface damage that might occur in the tank walls Thes e stress ^ 

particularly of concern during the initial cooldown period for a new tak p JcedTnTo 

sernoe, and dunng any subsequent thermal cycles. A preliminary thermal stress analysis 

of a high-pressure cryogenic storage tank was performed. Stn^ dSne TS 

operation were determined, as well as the trans.ent temperature distribution M elaSc 

analysts was used to determine the thermal stresses in the inner wa S Lff p'T 

temperature data. The results of this elastic analysis indicate that the toer wall of the 

s orage ta will experience thermal stresses of approximately 145 000 psi (1000 MPa) 

2. T", 1S ?». room-temperature yield strength of 304L “ 

■ ’ T T Z;NUUU pS1 (17 ° MPa >- For th 'S preliminary analysis several 

™304L^t£S T? T y " bee ” considerci These fetors include increiLd’ strength 
04L stainless steel at cryogenic temperatures, plastic material behavior and increased 

strength due to sham hardening. In order to more accurately ZL the “ 
shesses and them affect on the tank material, further investigation is required pLicuS v 
m the area of material properties and their relationship to stress. Pearly 

2. Introduction 

A newly constructed cryogenic storage tank will experience a significant thermal eradient 
dunng its initial cooldown. This process frequently uses liqutd dtrogen tfco a 
storage vessel used to store liquid oxygen (LOX), with the LOX beifg pumped In o the 
storage vessel after the tank is cooled sufficiently. Note that the saturation temperataL 

(180 R) fOT°LOX Th “ atm ° S ff C preSSUr= ^ 77 K C 1 39 R) &r LN2 and 100 K 
(180 R) for LOX. The effects of these thermal stresses on the condition of the tank 

sts. =■ =0£T‘ ” - •»« « - — ~ 

Tzszs: z ir, z. 

away as it encounters the hot inner wail of the tank. This combination of flowing bei ng 

Pr l af0 ff d - C0nVeCti011 boiIi "S con ^ibon on the surface of ie J? 
This condition is characterized by extremely large values of the heat transfer coefficient. 

Jff'TnfTl f 0pera, ‘ on . 0 f the tank consists of the tank being nearly filled with LOX at 
approximately atmospheric pressure. The typical ullage (amount of unfilled space in the 


storage tank) is 10%. Additional LOX must be supplied to the tank to replace the amount 

— r heat ieak into the system - ™ s st ° rage operation c ° m P risK *<= 

f i / 7 ■ ? 5 ' ogH “ < : storage tank life. Under these circumstances; the temperature 

o f the tank wall !S approx, mately uniform at the saturation temperature of the stored fluid- 
tor LUX at ambient pressure, the temperature is 100K (180 R). 

High-pressure operation of the tank occurs when the cryogenic fluid is being supplied for 
a . P articular task, rocket engine component testing for example. Under these 
circumstances, another gas, typically nitrogen, is introduced at high pressure into the 
ullage space at the top of the tank, thereby pressurizing the tank and its contents. 

This analysis will be performed for a tank whose diameter is 105 inches, wall thickness 

ot 13 inches and a maximum operating pressure of 8500 psi. The tank material is 304L 
stainless steel. 

3. Analysis 

The analysis of the stresses for the high-pressure cryogenic storage tank should include a 
e ermmation of the stresses under normal high-pressure operation at cryogenic 
temperatures, thermal stresses under the initial transient chilldown operation as the tank 

IS vn!f m \° t? 102 ’ f Wdl 83 &e effect of re P eated dermal cycles (warm up and 
chilldown). This analysis determines the spatial distribution of the operational stresses 

and a preliminary value of the thermal stress at the inner wall. Further analysis is 
required to more accurately model the thermal stresses and to estimate crack growth. 

3.1 Operational Stresses Analysis 

Although unrelated to the thermal stresses (since the cooldown and normal operation are 
separate and distinct events) an analysis of the stresses due to the pressurization of the 
ank was performed. The radial and tangential stresses were determined using equations 
or thick-walled vessels [1], These equations are for thick-walled spherical shells, with 
an imposed internal pressure, p, and have been modified for the nomenclature used hr this 
document. They do not mclude the effects of weight or support points. The tangential 
stresses (a 1 and ct 2) shown in Equation 1 are equal due to the symmetry' of the sphere. 
The radial stress (<t 3) is. determined using Equation 2, while the shear stress (t) is 
determined from Mohr’s circle as in Equation 3. The radial location within the spherical 
shell is given by r, while r* and r 0 represent the inner and outer radii respectively. 
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Eqn. 1 


Eqn. 2 


Eqn. 3 


A vessel is considered to be thick-walled if the ratio of the wall thickness to the inner 
radius is greater than about 0.1. For the tank described in this analysis, this ratio is about 
U.25, definitely a thick-walled vessel. For this analysis, a pressure of 8500 psi was used 
resulting in a maximum stress of -8500 psi (the negative sign indicates a compressive 



stress) in the radial direction at the inner wall of the tank. The maximum tangential stress 
was about 1 7,800 psi, while the maximum shearing stress was about 13,000 psi, both 
occurring at the inner walk These results indicate that the imposed pressure produces a 
maximum stress approximately equal to the 17,500 psi allowable stress specified by the 
ASME Code, Section VIII for. 304L stainless steel storage vessels at or below room 
temperature [2], The spatial stress distribution is shown in Figure 1. Note that the” 
tangential stresses are tensile stresses, while the radial stress is a compressive stress. 
Recall that these stresses are due to pressurization effects only; additional stresses due to 
the weight and support of the structure itself are not included. If desired these stresses 

could be determined and superimposed on the pressure stresses to find the total stress in 
the tank. 


It is interesting to compare these results to the yield and ultimate strengths of 304L 
stainless. ASTM sets a yield strength of 25,000 psi (170 MPa) and an ultimate strength 
of 70,000 psi (485 MPa) as minimum requirements for 304L stainless steel, while Yuri, 
et al • Pi have determined these strengths at LN2 saturation temperature (77 K / 139 R) to 
be about 50,000 psi (o50 MPa) and 225,000 (1550 MPa) psi respectively. 



Radius (in) 

Figure 1. Stress distribution in thick-walled spherical vessel (high-pressure operation). 

3.2 Thermal Analysis 

As the liquid nitrogen is pumped into the storage tank from the pipe at the bottom of the 
tank to cool the tank, the liquid nitrogen entering the tank will boil as it touches the inner 




walls, which are initially at ambient temperature. Because the incoming nitrogen is both 
owing and boiling, the heat transfer coefficient in this situation is quite high on the 
order of several thousand- it is this high- heat tiler coefficient 

rei ed r Ae r . e f VeIy J 0W £hermaI conductivity and the large spatial dimensions that 
result m the expected significant thermal stresses. 

A thm-walled tank, or a tank made from a high thermal conductivity material, would not 
experience such dramatic thermal stresses. The temperature would be much more 
uniform throughout the matenal at any given time during the cooldown procedure 
resulting m relatively little thermal stress, which is caused by the spatial thermal gradient.’ 

Temperature Distribution 


The temperature profile can be obtained by solving the heat conduction equation, a partial 
differential equation seen m Equation 4 for spherical coordinates. As can be seen this is 
a three-dimensional, transient problem. 
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Eqn. 4 


) th °. l i gh J the probIem 1S a three-dimensional transient conduction problem, it can be 
simplified to a one-dimensional transient conduction problem in the radial direction. The 
reasons for the simplification are the relatively low thermal conductivity value for 
stainless steel and the relative magnitude of the spatial dimensions of the tank. The 

•P r0 _ em , can , be fiirtiler simplified for this preliminary analysis by treating the 
tuCi mopnysicai properties, such as thermal conductivity, as constants. Although 
acceptable for this analysis, this is not physically accurate, as the thermal conductivity 
vanes with temperature over the range from ambient to cryogenic temperatures The 
simplified version of the heat conduction equation (ID, transient, constant properties) is 
shown below m Equation 5, where r is the radial location, T is the temperature, t is time 
and a is the thermal diffusivity. 


2 3T d*T 1 3T 

r 3r dr 1 ~~a~dt E P n - 5 


As with any partial differential equation, initial and boundary conditions are required to 
completely define the problem. Initially, the entire sphere is at ambient conditions, so a 
uniform temperature initial condition is appropriate, which is defined in Equation 6. The 
inner surrace of the sphere is subjected to forced-convection boiling, which is located 
only in the region occupied by the liquid, so the boundary condition for the inner wall is 
actually transient, depending on the fill rate or the amount of ullage in the tank. For this 
preliminary analysis, a steady boundary condition representing the forced-convection 
boiling is used. This is justified, since the area of interest in the initial part of the 
cooldown is located at the inlet, which is at the bottom of the tank. The portion of the 
tank not covered with LN2 has a much smaller heat transfer coefficient, and is being 
cooled much more slowly, generating much smaller thermal stresses. Therefore, as soon 
as there is liquid in the tank, the steady boundary condition is locally appropriate. The 
outer wall of the sphere is encased in a vacuum jacket, which contains insulation. For 
this preliminary analysis, this boundary will be modeled as an insulated surface. This is 



Temperature (R) 


acceptable for the modeling of the initial cooldown, but would need to be modified for a 
steady state thermal analysis of the tank to include the heat leak into the system. The 
eonveetioH—and- insutated— boundary conditions are shown - ^Equations 7 and” 8 
respectively, where t is time, r is the local radial location, and q and r 0 represent the inner 
and outer wall radii respectively. The heat transfer coefficient is given by h, and the 
thermal conductivity by k, both of which are found in Equation 7. 


T{r,t = 0) = 7) 

= h{T(r h t)-T„ L, 


dr 


dT_ 

dr 


= 0 , 


Eqn. 6 
Eqn. 7 


Eqn. 8 


The solution to this properly posed partial differential equation, which mathematically 
models the transient heat conduction in the sphere, can be approximated usins finite 
difference techniques. Specifically, central differences were used for the spatial 
derivatives, while the explicit method was used for the time derivative. The resulting 
spatial temperature distribution is shown for various values of time (starting with t=0 and 
ending with t=6000 seconds) in Figure 2. 



Figure 2. Temperature Distribution in Tank Wall for Selected Time Values. 




s ^ ow that the temperature gradient is quite steep at the beginning of the 
chi Idown process then tends to gradually decline in the first ten minutes ^er note 

cryogenic fithdln tST^T ^ eSSentiaiIy achieved the saturation temperature of the 
transfer coefficient T 

is°seen ITT °“ CT ^ 

j • . , r minutes (6000 seconds) the temperature of the outer wall is nnw 

ecreasing, avmg dropped approximately 50 degrees at this point in the process. 

Thermal Stresses 

^ the iim % Wa11 of the sphericaI vessel attains the saturation 
JSv r Uld (IiqUid nitr ° gen f0r the C001 down operation) 

^ inutes ’ J due P^Tto^i^rhS 

then^a. s^Ts flwr eqUa ‘'° n ' “ e daStic theimaI ^ to <hf 


*'=B (Ar) 


Eqn. 9 


coefficient ofTeS *%T° duIus of eIasti <%> v is Poisson’s ratio, and y is the 
the Z, i ° f , expansion - (^e subscript T denotes thermal stress.) The resulting 

rmal stress for this situation is approximately 145,000 psi, using the temperature 
difference for the inner wall of about 390 R (difference be^een LbLt ^N2 
I .° n te f 1 P er 2 hires) for a liquid nitrogen chilldown. This stress is due to the thermal 
f. nn f m Sphere ’ Where the lnner waI1 is at the LN2 saturation temperature while 

tenshe T 7h temperature. The thermal stress on the inner wall is 

tensile, as the cooling tends to contract the material. 

Typical values for the yield and ultimate strength of 304 L stainless step! ■>< nnn 

r 7 °’° 0< ! psi (485 MPa) ’ per ,he as™ ml. 25 bZZ 

„ ar ess resul '' " y ould seem tha t tie tank would be expected to fail However 

much taef (h»T ! Z ' ri eW and " itimate Stren S ths at ^genic temperatures are 
ucn larger than the room temperature values. So even as the temperature rn-*rhW 

creates thermal stresses, the cold-material has a much higher strength. P Also reSll that 

stLtes^ti^r and f° eS n0t inClUde the effects of plasticit yi the strength of 
stainless s tee at room temperatures increases dramatically due to strain hardening with 

maximum values (yield and ultimate) of 1 60,000 and 1 85,000 psi. 

3.3 Future Analyses 

The analysis results will depend on the effects of these factors, strength as a function of 

imer S^mSaF fsT? § ' A$ t0 previousl T’ if the temperature of the 

• M ! “ atenaI 1S 77 then there is a thermal stress of about 145,000 psi but the 
yield and ultimate strengths of 304L at this temperature are about 50,000 md 225,000 



uSm2 1 T^stL^ b r W0Uld n0t be »P-W to reach the 

sr- is — ■> « - — «: 

^idu" and strain ‘ wica,,y “>« a» 

thermal cooldown. is ^^fsZZfrTl^ "* m3tenal *“« d “™S *e 
thermal stresses will return to zero as there is no In aChleV f 3 ““ foml temperature, the 
However, if plastic deformation has occurred a nn!^ * to drive them - 

and into the sphere may be in a state of mm ’ ? 10n ° f * S matenaI at * e inner wall 

Should this be the casef this residual stZ iT ? S ? eSS due t0 the P iastic deformation, 
imposed by the internal pressure during °h‘ h gebraicaIly °PPOse the tensile stresses 

investigation of the material behavior is^equired^mTesnhiT supposition. ** 

4. Conclusions 

answer t^^ue^ioi^r^^din^^e^ateriafb^h made ’ investigation is required to 
and subsequent thermal cycles This ^ h ^1°^ unn S this initial cooldown process 
as well Z 3 search of the literature 

the effects of this initial cooldown u™ .dfL _® Imp0 , rt ™ t parameters and determine 


the effects of this initial cooldown. More fo„ , . 

as part ot this further mvestieartnn Ac *V;«. •* 7 J “* «« axso neeaea 

about the stress state of the material this info ' n '’. estl ® atl ° n Provides better information 
in the tank wall. ' S ,nfoimatron ca ” be nsed to model eraok growth 
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